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We calculate the phonon density of states (DOS) for strongly amorphous materials with a short- 
ranged interatomic potential. Exponentially decaying and abruptly truncated interatomic potentials 
are examined. Thermally excited mean square deviations from equilibrium are calculated with rapid 
increases noted as the average number of neighbors is reduced. The Inverse Participation Ratio (IPR) 
is used to characterize the phonon states and identify localized phonon modes as the bonding range 
(and hence the average number of neighbors per atom) is diminished. For the truncated potential, 
the characteristics of the IPR histogram change qualitatively below nncigh with the appearance of 
localized phonon modes below nncigh = 6.0. 

PACS numbers: 63.20.-c, 63.20.Pw, 63.50.Lm, 63.50.Gh 



INTRODUCTION AND METHODS OF 
CALCULATION 



The range of materials used in current engineering ap- 
plications is diverse, and many structural media have 
an amorphous character. Acoustic and thermodynamic 
properties relevant to the performance of a class of mate- 
rials are technologically significant properties. Wc exam- 
ine highly amorphous systems and calculate the phonon 
frequency density of states (DOS) which may be used to 
determine the thermally excited deviations from equilib- 
rium. In addition, the Inverse Participation Ratio (IPR) , 
i^YPK — provides a way to characterize 

the phonon modes and determine the degree to which 
they are localizedii^. In the moment ratio kjpr, ip is the 
positional amplitude of a phonon state, and vibrational 
modes with localized character are associated with a fi- 
nite value of KipR whereas for extended states kipr tends 
to zero in the bulk limit. 

An important question is the effect of changing the 
range of the coupling between atomic species. A short- 
ranged coupling corresponding to a relatively loosely 
packed system will lead to reduced rigidity, while a longer 
ranged coupling increases the average number of nearest 
neighbors, and thereby imparts rigidity to the amorphous 
lattice. 

We consider a coupling scheme with an abrupt cutoff 
Tc, and we also examine an exponentially decaying in- 
teraction scaling as e~^'^ where the rapidity of the decay 
may be tuned with a finite (though rapidly diminishing) 
coupling as a function of distance by changing the damp- 
ing constant A. Hence, with genuine discontinuities ab- 
sent from the interaction function, qualitatively similar 
(albeit less abrupt) changes will occur as the parame- 
ter controlling the rapidity of decay is modified. The 
flat potential terminating beyond a threshold radius 
is examined in Section II, and in Section III results are 
discussed for the exponential coupling scheme. 

In analyzing lattice dynamics, we use the classical La- 



grangian formalism 

i—l 2—1 j — 1 

where the factor of 1/2 in the potential is included to 
compensate for redundant counting of bond energies, rii 
specifies the number of neighbors for the atom given the 
label z, and Vij(r) is the specific interaction between 
atoms i and j. In the context of our calculation, the 
system is contained in a supercell with cubic geometry 
and length L, containing TV atoms of equal mass m with 
periodic boundary conditions assumed. 

In the temperature regimes of interest, we assume the 
deviations in atomic positions to be small in relation to 
the equilibrium bond lengths , and we employ the har- 
monic approximation (i.e. essentially a Taylor expansion 
of a well behaved potential) with 
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with Kij being the second derivative of V{rij) at 
In spite of the implementation of the harmonic approx- 
imation in the inter-atomic bonding potential, the non- 
coUinearity of bonds will still give rise to anharmonic 
terms in the potential. 

Nevertheless, a further harmonic approximation may 
be appropriate if the atomic coordinates are written in 
terms of the equilibrium coordinate with a shift term, 
such + Sf with similar expressions for the y 

and z coordinates; the bond length lij has the form 
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and A°f 

''J 



(z? — z°). Hence, the potential energy stored in the 
strongly disordered lattice depends only on the difference 
of coordinates such as, e.g., for the equilibrium x co- 
ordinates and {Sf — Sj) for the corresponding shifts from 
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equilibrium. If the latter are sufficiently small relative 
to the former, it is appropriate to expand about A^'^, 

A^/, and AO/, yielding (k, ~~ l%f « [a,, • {1-5,) 

where A,;j is a unit vector formed from the difference 
A^j = X? — x°. The interatomic potential, to quadratic 
order in components of 5i and 5j is 
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i=l 3 = 1 



The standard Lagrangian formalism then yields the equa- 
tions of motion, given by, e.g. + = 0. Written 

in terms of the vector displacement J^, one has (we as- 
sume the atoms to have identical mass to) 
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Using the ansatz bi = 5i e*"* removes the explicit time 
dependence and reduces the solution of the equations of 
motion to an eigenvalue problem, which may be written 
as 
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The phonon modes are obtained by diagonalizing the 
symmetric matrix, and in order to be definite, the eigen- 
states are taken to be normalized to unity. The form of 
the matrix is illustrated in Fig. [1] for a system contain- 
ing four particles and rendered schematically below the 
matrix on the far left. Immediately to the right are neigh- 
bor connectivity graphs illustrating the bonding scheme 
which gives rise to the matrix. The latter provide a way 
to set up the nonzero 3x3 submatrices away from the 
main diagonal. 

In calculating thermodynamic quantities such as the 
mean square deviation from equilibrium per atom JrmSj 
one evaluates the partition function 
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(e-'^")dconf. 
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where /3 = 1/A:b with fee being the Boltzmann constant; 
the integral sign and "dconf" refer to summing over all 
possible system configurations, where all possible atomic 
and velocities deviations are sampled. One may deter- 
mine the mean square deviations from equilibrium with 
the aid of formalism developed elsewhere^. The RMS 
displacements from equilibrium (5rms j &t the level of the 
harmonic approximations we have made will consist of a 



FIG. 1: (Color online) A sample matrix is shown for a four 
member system depicted schematically below the matrix. The 
connectivity graphs illustrate which neighbors correspond to 
each of the numbered particles, and the connectivity pattern 
dictates the form of the matrix to be diagonalized to obtain 
the phonon modes. 



thermal factor proportional to T^/^ and a factor depen- 
dent on the charateristics of the lattice. The latter be of 
most interest in this context, because it will depend on 
the rigidity associated with the specific atomic configu- 
rations and the bonding pattern among the atomns com- 
prising the amorphous material. This normalized mean 
square deviation lacks the thermal factor, and is propor- 
tional to the square root of the sum over the reciprocals 
of the squares of the phonon frequencies (with the zero 
frequencies excluded). Hence, the normalized RMS dis- 
placements i5J^f^^g are given by 



-•RMS 




(8) 



where M is the total number of nonzero frequencies. As 
the system becomes less rigid and zero modes prolifer- 
ate, the lattice will lose local stability, ultimately being 
reduced to small isolated clusters which may move inde- 
pendently. As the integrity of the amorphous medium is 
lost, the notion of RMS displacements from positions in 
a single large lattice ceases to be an appropriate quantity 
to consider. 

For a specific configuration of disorder, one calculates 
the phonon density of states by diagonalizing the appro- 
priate matrix to determine the frequencies of the 3N 
phonon modes. To build up a histogram with sufficient 
statistics, it is important to sample multiple realizations 
of disorder, and the total number of eigenvalues used in 
constructing histograms will be 3(iV)nconf where {N) is 
the mean number of atoms in the unit cell, and nconf 
is the total number of configurations sampled. To en- 
sure that at least phonon frequencies are sampled. 
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we set riconf = ^{n-u;/ i^))- In calculating the DOS his- 
tograms we have n^^ = 10^, whereas for the IPR calcula- 
tions where both the phonon eigen-frequencies and states 
must be obtained, we sample = 10^ eigenstates. 

In the highly amorphous arrangements of atoms we 
consider, atomic positions are uncorrelated, and hence 
there is no local order. The probability of having n 
atomic members in a volume v = is p{n) = ^^^i^ e"^" 
with p being the number density. We sample realiza- 
tions of disorder in an unbiased way with a stochastic 
technique based on the Metropolis Criterion. Beginning 
with the integer closest to the mean occupancy (n) = pv, 
we make a scries of attempts to change n (the number 
of attempts is made equal to n to ensure proper rando- 
minization) where half of the attempts would reduce n 
by one, and half try to increment n by one unit. For 
attempts to increase n, the relevant parameter is the 
probability ratio r+ = p{n + l)/p{n) = pv/(n + 1), and 
the increase is accepted if X,- < r+ where X,- is a ran- 
dom number confined to the interval [0, 1] and generated 
with uniform probability. On the other hand, the ap- 
propriate probability ratio for attempts to decrease n is 
r_ =p{n—l)/p(n) = n/pw, and the occupancy is reduced 
by one if X^ < with X,- a random variable again sam- 
pled uniformly from the interval [0, 1]. In generating the 
random numbers, a Mersenne twister algorithm is used 
to minimize correlations effects between successively gen- 
erate numbers and to ensure a long period in the random 
number sequence. 

With the occupancy n determined, equilibrium coor- 
dinates y°, and z° for each atom are chosen with 
uniform probability in the interval [0, L], and the appro- 
priate Hamiltonian matrix is constructed and diagonal- 
ized to calculate the photon frequencies to construct the 
vibrational DOS; the eigenstates themselves are retained 
to calculate the Inverse Participation Ratio (IPR). To re- 
duce the impact of finite size effects, we assume periodic 
boundary conditions. 

To obtain good statistics and reach the bulk limit, 
we calculate the phonon frequency spectra for multiple 
disorder realizations for several different system sizes. 
Averaging over disorder provides a smoother density of 
states curve in a manner which approximately mimics the 
DOS profile in the thermodynamic limit where some self- 
averaging would be expected to occur. In the sequence 
of successively larger systems examined, the merging of 
the curves indicate the attainment of the bulk limit. 

To illustrate the atomic configurations and lattice con- 
nectivity for the highly amorphous material with a trun- 
cated coupling between atoms, snapshots of a portion 
of the strongly disordered medium are shown in Fig. [2l 
Fig. [3J and Fig. |31 Although the atomic configurations 
are the same in each case, the different values of the trun- 
cation radius Tc yield very different bonding configura- 
tions. The first image in Fig. [2] shows the medium with 
'T-ncigh = 3.0 and hence a very low connectivity among the 
atoms. Fig. [3] is an intermediate case where rinoigh = 6.0. 
The last of the depictions in Fig. 0] is a configuration 




FIG. 2: (Color online) Plot of atomic locations and bond- 
ing scheme for rincigh = 3.0. Red symbols represent atoms, 
and the dark lines connecting the spheres are bonds between 
neighbors. 




FIG. 3: (Color online) Plot of atomic locations and bond- 
ing scheme for ^ncigh = 6.0. Red symbols represent atoms, 
and the dark lines connecting the spheres are bonds between 
neighbors. 



where the average number of neighbors is relatively high, 
and rineigh = 12.0. The bonding scheme and its effective- 
ness in setting up a rigid lattice varies significantly with 
increasing nnoigh- 



II. TRUNCATED INTERACTION WITH A 
CUTOFF 

We calculate the phonon density of states curves for 
a range of values of Tc- As the threshold radius is de- 
creased, the average number of neighbors rinoigh for each 
atom also is decreased, and it is less likely the amor- 
phous configurations will be rigid; the migration of the 
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FIG. 4: (Color online) Plot of atomic locations and bond- 
ing scheme for rineigh = 12.0. Red symbols represent atoms, 
and the dark lines connecting the spheres are bonds between 
neighbors. 
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FIG. 5: (Color online) Normalized phonon Density of States. 
Panel (a) corresponds to Tincigh = 10.0, panel (b) displays 
DOS curves for rinoigh = 15.0, and panel (c) is plotted for 
»^neigh = 20.0. Results for various supercell size L are shown. 
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FIG. 6: (Color online) Residual histogram amplitudes for 
iT-ncigh = 10.0, Tincigh = 15.0, Tincigh = 20.0 where uumerals 
enclosed in boxes indicate the systems compared; e.g. "12- 
14" refers to the residual amplitude between systems with 
i = 12 and L = 14. 
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FIG. 7: (Color online) Normalized phonon Density of States. 
Panel (a) corresponds to Tinoigh = 10.0, panel (b) displays 
DOS curves for Tindgh ~ 15.0, and panel (c) is plotted for 
TT'ncigh = 20.0. Results for various supercell size L are shown. 



main peak of the roughly uni-modal phonon density of 
states toward lower frequencies is compatible with this 
expectation. However, one must be sure that L is large 
enough that finite size effects have minimal impact on 
the DOS curves. Figure [5] shows the density of states 
curves for a range of system sizes, and the near coinci- 
dence of the DOS profiles for the larger systems is ev- 
idence of convergence with respect to the size L of the 
supercell. Figure |6] highlights the difference among the 
DOS curves and depicts the residual differences between 
the disorder-averaged curves corresponding to different 
values of L. The residual values, differences between DOS 



curves for successive supercell size L shown in panels (a) , 
(b). and (c) of Fig. [S] are two orders of magnitude below 
the principle phonon density of states amplitude, pro- 
viding additional evidence the DOS curves arc converged 
with respect to the supercell size L. 

The principle results are summarized in Fig. rcf- 
fig:Fig7, where the DOS curves for a range of average 
neighbors (i.e. from rineigh = 10.0 to rineigh = 20.0) are 
graphed for the largest system size (L = 14.0) in the 
study. As the average number of neighbors is decreased 
from rinoigh — 20.0, the shape of the curve changes and 
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FIG. 8: (Color online) Fraction of zero frequency phonon 
states for system sizes ranging from L = 4.0 to L = 10.0 
with the mean number of neighbors, rinoigh on the horizontal 
axis. 
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FIG. 9: (Color online) Normalized mean square deviations 
plotted versus Tincigh for various supercell sizes L. Panel (a) 
is a closer view of the rapid jump in the RMS fluctuations 
near rinoigh = 10.0, and panel (c) show the mean square fluc- 
tuations with respect to L for several relatively large values 
of 

^neigh . 



the weight of the histogram shifts leftward, toward the 
lower frequency regime. Ultimately, in the vicinity of 
nneigh = 10.0, a sharp zero frequency peak forms on the 
left edge of the graph, and increases in height with de- 
creasing Tc. 

Fig. [5] shows the relative fraction / of zero modes with 
respect to the mean number of nearest neighbors rincigii. 
There is an abrupt increase in / beginning in the vicin- 
ity of rinoigh = 8.0 with an approximately linear increase 
in the zero mode fraction with decreasing rinoigh begin- 
ning for rinoi-^h < 8.0. The rapid proliferation of zero 



modes is a phenomenon which one may attribute to the 
degradation of rigidity with decreasing average number 
of nearest neighbors rincigh- Fig. [HI shows the normalized 
(the t^^"^ factor is not included) root mean square (RMS) 
deviation which would be calculated thermodynamically 
as discussed in Section I. Again, an abrupt (albeit noisy) 
increase in (S^I^^g begins in the vicinity of rinoigh = 8.0. In- 
set (a) of Fig. [9] is a closer view of the 5^^^^ curves near 
the rapid increase. As the average number of neighbors 
'T-noigh is decrease, the mean square deviation continues to 
increase, ultimately reaching a maximum in the vicinity 
of 

^ncigh — 5.0. The rapid decline which begins in i5J^]y[g 
for nnoigh > 5.0 is not an indication that the amorphous 
lattice has become more rigid, but is instead a result of 
the rapid proliferation of zero modes which indicate the 
rigidity of the medium has been compromised and the 
RMS deviation from equilibrium is no longer a reliable 
index of the integrity of the lattice. 

Inset (b) of Fig. [9] show the smoothly convergent be- 
havior of the thermally excited RMS deviations for rela- 
tively large values of the average number of atomic neigh- 
bors where there are very few zero frequency phonon 
modes. In this regime of relatively high rincigh values, 
fluctuations about equilibrium positions in the lattice are 
ultimately stable with increasing L and the specific con- 
figuration of atomic positions (albeit strongly disordered) 
is stable at finite temperatures in the bulk limit. 

The inverse participation ratio (IPR) is a useful index 
of the characteristics of phonons with respect to localiza- 
tion. The diagnostic merit of kipr is realized in the bulk 
limit where the IPR either converges to a finite value (for 
localized states) or tends to zero (for extended states). 
We investigate the extent to which phonons are localized 
by calculating the probability density values and aver- 
aging over a sufficient number of disorder configurations 
(i.e. we consider at least 10^ phonon states). Results are 
shown in Fig. [101 HH f^d [12] for supercell sizes L ranging 
from L = 4.0 to L = 10.0. 

The three sets of figures each contain four IPR his- 
togram graphs, where rinoigh is increased from one panel 
to the next, and between successive figures. In partic- 
ular, rinoigh ranges from rineigh = 2.0 to rineigh = 5.0 in 
Fig.[Tni from rineigh = 6.0 to rinoigh = 9.0 in Fig. [TTl and 
from rinoigh = 15.0 to rinoigh = 20.0 in Fig. [T^j where a 
relatively large average number of neighbors are involved. 

The graphs in Fig. [TUl Fig. [TTl and Fig. [H illustrate 
salient trends in the characteristics of the phonon states 
with respect to localization. IPR curves in panels (a), 
(b) , (c) , and (d) of Figure [TOl have a bimodal structure 
where there is a peak for relatively high IPR values as 
well as a sharp maximum at lower values of the Inverse 
Participation Ratio. The two peaks behave very differ- 
ently as the system size L is increased; the portion of the 
histogram near the rightmost peak values tend to con- 
verge, ultimately shifting neither in position, weight, or 
overall form with increasing L. On the other hand, the 
peak appearing at lower IPR values continues to sharpen 
and migrate leftward, toward even lower IPR values as 



6 



mplitudt 

CO 

b 






< 




B2.0 








o 








N 




B 

go.o 






IL = 12.01 



IL = lO.Ol 



"neigh = 3.0 



pM] (b) 

^ IL =6.0 1 

=4.0| 




■§5.0 

I4.O 

13.0 

i2.o 

1 1-0 
go.o 





L =4.0 — 
L =6.0 — 
L=8.0 — 
n L = 10.0 — 




It 



•S4.5 





4.0 







L =4.0 — 
L =6.0 — 
L =8.0 — 
Z. = 10.0 — 
A L = 12.0 — 
M t=14.0 — 


\ n . ^=20.0 

1 neigh 









-3.0 -2.0 

Log,„(IPR) 



-1.0 



0.0 



FIG. 10: (Color online) IPR histograms for various system 
sizes with rinoigh = 2.0 in panel (a), nncigh ~ 3.0 in panel (b), 
?^neigh = 4.0 in panel (c), and Jincigh = 5.0 in panel (d). 



FIG. 12: (Color online) IPR histograms for various system 
sizes with n,icigh = 15.0 in panel (a), rineigh = 17.0 in panel 
(b), Tincigh = 19.0 in panel (c), and rineigh = 20.0 in panel (d). 




FIG. 11: (Color online) IPR histograms for a range system 
sizes with rineigh = 6.0 in panel (a), rineigh ~ 7.0 in panel (b), 
rineigh = 8.0 iu panel (c), and rineigh = 8.0 in panel (d). 

L is made larger. The steady shift toward lower Inverse 
Participation Ratios is particularly significant given the 
fact that the horizontal axis in the graph is a base ten 
logarithm of the IPR; the substantial shift of density to 
lower IPR values with the IPR quickly shrinking signals 
the presence of extended states where finite size effects 
interfere less with the extended character of the corre- 
sponding phonon modes. 

For the right-most peak appearing in the higher IPR 
regime behaves in a qualitatively different way with in- 
creasing L, eventually ceasing to move or change in shape 



with further enlargement of the supercell. The latter in- 
dicates the presence of localized states where the phonon 
modes are characterized by a finite IPR value. The con- 
vergence of the histogram curves with respect to system 
size may be attributed to the diminution of the severity 
of finite size effects as the regime where L ^ ^ is at- 
tained, with ^ being the localization scale of the phonon 
modes. As rineigh is increased within the figure, the bi- 
modal structure remains, but the convergence of the peak 
in the large IPR regime is less rapid and most of the sta- 
tistical weight eventually shifts toward smaller IPR val- 
ues, or in the direction of extended states. 

In Fig. [TT] there is a sudden change in the characteris- 
tics of the Inverse Participation Ratio curves where the 
bimodal profile gives way to a single peak structure. The 
transition occurs for approximately the same value of 

^neigh 

where the RMS deviation (5rms peaks in Fig. [HI 
In all cases shown, with rineigh ranging from rineigh = 6.0 
to rineigh = 9.0, thc peak becomes sharper and migrates 
toward lower values of the Inverse Participation Ratio. 

Finally, results for the regime of large n„eigh are shown 
in Fig. [T21 In particular, panels (a), (b), (c), and (d) 
of Fig. [12] correspond to average neighbor values ranging 
from rinoigh = 15.0 to rinoigh = 20.0. A salient charac- 
teristic of the IPR histograms is the lack of significant 
variation and location of the peaks for a fixed supercell 
size L. 

The apparent approach of the Inverse Participation 
Ratio histogram curves to limiting profiles suggests that 
a strongly rigid limit has been achieved for rinoigh ^ 15.0. 
Statistically speaking, in this Tc limit, lattice vibration 
modes have identical characteristics with respect to lo- 
calization. 
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FIG. 13: (Color Online) Normalized phonon Density of 
States. Panel (a) corresponds to A = 2.0, panel (b) displays 
DOS curves for A — 1.0, and panel (c) is plotted for A — 0.5. 
Results for various supercell size L are shown. 




FIG. 14: (Color Online) NormaUzed phonon Density of 
States. Panel (a) corresponds to A = 2.0, panel (b) displays 
DOS curves for lambda — 1.0, and panel (c) is plotted for 
A = 0.5. Results for various supercell size L are shown. 



III. EXPONENTIALLY DECAYING 
INTERACTION 

An exponentially decaying interaction scheme provides 
a short-ranged interaction without the need for an abrupt 
cutoff Tc between atoms is taken to be proportional to 
e~^^ where 1/A provides a length scale. The DOS curves 
shown in Fig. [T3] are obtained for a range of decay con- 
stants A. Panels (a), (b), and (c) of Fig. [T^ correspond 
to decay constants A = 2.0, A = 1.0, and A = 0.50. In 
seeking the bulk limit, two relevant length scales are the 
typical separation I = p~^l^ (unity in the present case 
where p = 1) between atoms and A^^ for the interatomic 
potential. Hence, one must have L ^ max A~^ } 

for convergence to the bulk limit. 

The relatively rapid convergence of the phouou DOS 
curves for A = 2.0 is compatible with the short range of 
the interaction between atoms. For the relatively long- 
range cases A = 1.0 and A = 0.50, convergence is less 
rapid with the slowest approach to the bulk limit occur- 
ring in panel (c) of Fig. 1131 and a slightly more rapid 
attainment of the thermodynamic limit in the intermedi- 
ate case displayed in panel (b). 

As in the case of the interaction with a cutoff scale r^,, 
we calculate IPR histograms for the exponentially de- 
caying inter-atomic coupling, with results appearing in 
Fig. [M] The changes in the characteristics of the states 
with respect to localization are less dramatic with respect 
to localization with changes in A^^ than with different 
values for the case of the truncated potential. In pan- 
els (a), (b), (c), and (d) of Fig. [Ml the IPR histogram 
is shown for coupling decay constants A = 1.0, A = 1.5, 
A = 2.0, and A = 2.5. If the interaction is relatively 
long-ranged (i.e. comparable to or greater than the inter- 



atomic separation I = p^^^'^), very little change occurs in 
the shapes of the histogram curves for a particular system 
size L. The IPR histogram profiles are uni-modal, and 
the tendency of the peaks to become narrower and shift 
in the direction of lower IPR values is similar to what is 
seen for the IPR histograms corresponding to the trun- 
cated potential in the limit that the average number of 
neighbors is large. With shorter ranged and more rapidly 
decaying couplings between atoms, changes in the IPR 
histograms are still relatively subtle, in contrast to the 
abrupt shifts occurring in the case of the sharply trun- 
cated potential as Tc is reduced. A feature which does ap- 
pear in the case of the more rapidly decaying potentials 
[e.g. for A = 2.0 shown in panel (c) of Fig. [HI and A = 2.5 
in panel (d) of Fig. [Tl] , is a small shoulder on the right 
side of the graph for relatively large values of the Inverse 
Participation Ratio. Although the structure is relatively 
weak, it persists at approximately the same amplitude 
with increasing L, suggesting some of the phonon states 
are at least quasi- localized if the length scale A~^ of the 
exponential coupling is sufficiently small. This behavior 
is qualitatively similar, though markedly subtler, than 
the convergence of the right-most peak corresponding to 
localized states for the case of the truncated potential 
where nnoigh < 5.0. 

IV. CONCLUSIONS 

We have examined characteristics of phonon states in 
strongly disordered media where the coupling between 
atoms is finite in range, with either a fiat profile ter- 
minating at a radius fc, or a more gradually decaying 
exponential coupling scheme with its own length scale 
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A~^. In the case of the abruptly truncated interaction, 
an interatomic potential profile that may be more ap- 
propriate for relatively loosely packed particles with a 
non-covalent more rapidly decaying Lennard- Jones type 
of coupling, there are sudden transitions in shapes of the 
phonon DOS curves and the root mean square deviations 
<5rms from equilibrium with decreasing re- In addition, 
the Inverse Participation Ratios change qualitatively be- 
tween nnoigh = 5.0 and rinoigh = 6.0, converting from a 
bimodal profile with localization signatures to a single 
peak which migrates to lower IPR values with increasing 
supercell size as expected for extended states. 



For the exponential coupling profile, which decreases 
gradually with the inter-atomic distance, the variation 
of the IPR histogram with decreasing range is not 
punctuated by sharp changes. Moreover, significant frac- 
tion of the phonon states have extended character even 
for the relatively short-ranged case A = 2.5. 
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